Resonance phenomenon for the Galerkin-truncated Burgers and Euler equations 
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It is shown that the solutions of inviscid hydrodynamical equations with suppression of all spatial Fourier 
modes having wavenumbers in excess of a threshold K G exhibit unexpected features. The study is carried 
out for both the one-dimensional Burgers equation and the two-dimensional incompressible Euler equation. For 
large Kg and smooth initial conditions, the first symptom of truncation, a localized short-wavelength oscillation 
which we call a "tyger", is caused by a resonant interaction between fluid particle motion and truncation waves 
generated by small-scale features (shocks, layers with strong vorticity gradients, etc). These tygers appear when 
complex-space singularities come within one Galerkin wavelength A G = 2ir/K G from the real domain and 
typically arise far away from preexisting small-scale structures at locations whose velocities match that of such 
structures. Tygers are weak and strongly localized at first — in the Burgers case at the time of appearance of 
the first shock their amplitudes and widths are proportional to K a ~ 2 ^ 3 and K G ~ X ^ 3 respectively — but grow 
and eventually invade the whole flow. They are thus the first manifestations of the thermalization predicted by 
T.D.Lee in 1952. 

The sudden dissipative anomaly — the presence of a finite dissipation in the limit of vanishing viscosity after 
a finite time t*— , which is well known for the Burgers equation and sometimes conjectured for the 3D Euler 
equation, has as counterpart, in the truncated case, the ability of tygers to store a finite amount of energy in the 
limit A'g — > oo. This leads to Reynolds stresses acting on scales larger than the Galerkin wavelength and thus 
prevents the flow from converging to the inviscid-limit solution. There are indications that it may eventually be 
possible to purge the tygers and thereby to recover the correct inviscid-limit behaviour. 

PACS numbers: 05.20.Jj, 05.45.-a, 47.27Jv 



I. INTRODUCTION AND FORMULATION 

When the motion of a fluid is described at the microscopic 
level, a conservative Hamiltonian formulation is appropriate 
and statistical steady states can be described using Gibbs en- 
sembles. At the macroscopic level, however, one obtains a 
dissipative hydrodynamical description because macroscopic 
motion can be irreversibly degraded into thermal molecular 
motion. Curiously, T.D. Lee observed that, starting from the 
hydrodynamical or magnetohydrodynamical equations for an 
ideal fluid, one can obtain a conservative dynamical system 
to which Gibbsian statistical mechanics becomes applicable 
Ut]. For this he used a Galerkin truncation of the equations, 
a procedure that keeps only a finite number of spatial Fourier 
harmonics. For the case of the Galerkin-truncated 3D incom- 
pressible Euler equation, Lee obtained thermalized equilib- 
rium statistical states having an equipartition of kinetic energy 
among all the Fourier harmonics and thus a k 2 energy spec- 
trum. This is very far from the spectrum of fully developed 
turbulence as observed experimentally, which could lead one 
to believe that Galerkin truncation applied to the Euler equa- 
tion cannot tell us anything about the dissipative states of tur- 
bulence. 

Kraichnan was the first to think otherwise. Considerations 
of the Galerkin-truncated equilibria of the 2D Euler equation 
played an important role in his conjecture about an inverse 
energy cascade 10, [H. In 1989 he and S. Chen went much 
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further and wrote (fl, p. 162): 

the truncated Euler system can imitate NS 
[Navier-Stokes] fluid: the high-wavenumber de- 
grees of freedom act like a thermal sink into 
which the energy of low- wave-number modes ex- 
cited above equilibrium is dissipated. In the limit 
where the sink wavenumbers are very large com- 
pared with the anomalously excited wavenum- 
bers, this dynamical damping acts precisely like 
a molecular viscosity. 

Supporting evidence was found in 2005 with very-high- 
resolution spectral simulations of the 3D Galerkin-truncated 
Euler equation that showed the following: when initial con- 
ditions are used that have mostly low-wavenumber modes, 
the solutions have long-lasting transients in which only the 
high-wavenumber modes are thermalized, while the lower- 
wavenumber modes behave in a way similar to that for viscous 
high-Reynolds-number flow 15J. This seems to hold not only 
when the low-wavenumber modes are weak (as implicitly as- 
sumed by Kraichnan who invoked the fluctuation-dissipation 
relations) but also in the strong turbulence regime that dis- 
plays a K41-type inertial range. One possible interpretation is 
that the thermalized modes act as a kind of artificial molecular 
world, thereby allowing dissipative (Navier-Stokes) dynamics 
for the lower-wavenumber modes. 

We understand far too little about the mathematics of the 3D 
Euler and Navier-Stokes equations to start a serious analyti- 
cal investigation of what happens to solutions of the Galerkin- 
truncated 3D Euler equation when K G — > oo. However, such 
matters may be within reach for the one-dimensional invis- 
cid Burgers equation, a well-understood problem in the ab- 
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sence of truncation. Even in that "simple" case, the behavior 
at large K G is far from obvious. Indeed, there are known in- 
stances where an energy-conserving modification of the invis- 
cid Burgers equation with a small parameter is found not to 
converge l35ll to the inviscid limit Ja,01 (see also 01]). Hence 
caution is needed and we shall discuss this issue further. 

There is also an important practical reason to be inter- 
ested in Galerkin-truncated hydrodynamics. Spectral methods 
(and their pseudo-spectral variants) are among the most pre- 
cise methods for the numerical integration of hydrodynamical 
equations |9fl. By necessity, a finite resolution must then be 
used. In other words one integrates not the full hydrodynam- 
ical equations but their Galerkin-truncated modifications. If 
the high-wavenumber modes are sufficiently damped by vis- 
cous dissipation the difference may be extremely small. Yet 
the desire to push the Reynolds number can lead to serious 
truncation errors. Furthermore in investigations of the blow- 
up problem for the 3D Euler equation (cf., e.g., Refs. ITolfTTll 
and references therein) it is important to be able to distinguish 
genuine blow-up from truncation effects. 

We had one additional reason to investigate what exactly 
are the consequences of Galerkin truncation. M.E. Brachet 
(private communication 2007) informed us about a strange 
phenomenon observed when Galerkin truncation is used in 
conjunction with the one-dimensional inviscid Burgers equa- 
tion 

dfU + ud x u = 0; u(x, 0) = uq(x). (1) 

The initial condition uq (x) has just a few Fourier harmonics 
and the number of retained Fourier modes is large. The first 
symptom of Galerkin truncation found by Brachet was a spu- 
rious oscillation in physical space, seemingly born not at all 
where one would expect it, namely in the neighborhood of 
genuine small-scale structures such as shocks and their pre- 
cursors called preshocks, but completely "out of the blue" at 
a place having no particular small-scale activity, as illustrated 
in Fig.Q] 

In this paper we shall understand why this happens, both for 
the Burgers equation and for the Euler equation (so far mostly 
in 2D). The phenomenon is here called tyger after William 
Blake's poem for reasons given later. Before proceeding to 
explain the organization of the paper, it is useful to define 
our Galerkin-truncated problem more precisely for the case 
of the Burgers equation (the 2D Euler case will be formulated 
in Sec. Ill C\ . We restrict ourselves to 2ir periodic solutions 
which can be expanded in a Fourier series 

u(x)= J2 cikX{lk - ( 2 ) 

fc=0,±l,±2... 

Let K G be a positive integer, here called the Galerkin trunca- 
tion wavenumber. We define the Galerkin projector P^ g as the 
low-pass filter which sets to zero all Fourier components with 
wavenumbers \k\ > K G . In other words 

P KG u(x)= ]T c ikx u k . (3) 

\k\<K G 



The (untruncated) inviscid Burgers equation, written in con- 
servation form, is 

d t u + d x (u 2 /2) = 0; u(x,0)=uo(aO. (4) 

The associated Galerkin-truncated (inviscid) Burgers equation 
whose solution is denoted v(x, t) is obtained by applying the 
low-pass filter to both the initial condition and to the nonlinear 
term fl5ll . It reads 

d t v + P KG d m (v 2 /2) = 0; v = P KG u . (5) 

As is well known, the gradient d x u of the solution to the invis- 
cid Burgers equation with smooth initial data typically blows 
up after a finite time . At the solution u has a cubic-root 
singularity, called a preshock lfT6llT7ll . Beyond i* the solution 
can be continued by introducing a small viscous term into the 
r.h.s. of ([T); in the limit of vanishing viscosity one obtains 
what we here call the inviscid-limit solution, which has one 
or several shocks lfl2ll . This is a generalized solution which 
satisfies the inviscid Burgers equation only in a weak sense. 
The inviscid-limit solution has a dissipative anomaly, i.e., it 
dissipates energy even in the limit of vanishing viscosity. In 
contrast, the solution to the Galerkin-truncated equation (O 
stays smooth and conserves energy forever. 

This paper has two main parts: Sec. [II] deals with the nu- 
merical exploration of the tyger phenomenon and includes 
soft phenomenological interpretations of our various find- 
ings. More specifically, in Sec. Ill Al we identify the resonant 
particle-wave interaction mechanism responsible for the birth 
of tygers. In Sec. lIIBI we present the whole temporal scenario 
from the birth of tygers to full thermalization. In Sec III CI we 
show that the 2D incompressible Euler equation also gives 
rise to tygers by a mechanism similar to what we find for the 
Burgers equation. In Sec. Ill Dl we investigate the energetics 
(dissipative anomaly) and the issue of the (weak) limit of the 
truncated solutions when K G — > oo. Sec. [ITTlis restricted to 
the birth of tygers for the Burgers equation; it involves state- 
of-the-art simulations of the scaling properties with K G up to 
40, 000 and also a fair amount of analytic theory. Open prob- 
lems and conclusions are presented in Sec. [TV] There are four 
technical appendices. 

II. SIMULATIONS AND PHENOMENOLOGY 

A. Tygers and resonance 

Henceforth, when we write about a/the "untruncated solu- 
tion of the Burgers equation", without specifying more, it is 
the inviscid limit of the untrucated Burgers equation which is 
understood. A particularly simple 27r-periodic solution of the 
Burgers equation is obtained with the initial condition 

uq(x) = sin a;, (6) 

which has two stagnation (zero-velocity) points, x = with 
positive strain (gradient) and x = it with negative strain. The 
latter gives rise to a cubic -root preshock singularity at the time 
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t = i* = 1 Bal . Hereafter the initial condition (O and 
its space-translates will be referred to as "single-mode initial 
condition". 
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FIG. 1: (Color online) Growth of a tyger in the solution of the invis- 
cid Burgers equation with initial condition vq(x) — sin(x — n/2) 
(to avoid graphical edge effects). Galerkin truncation at K G — 700. 
Number of collocation points N — 16, 384. Output time as labelled. 
Observe that the bulge appears far from the place of birth of the 
shock. In the PDF online version of this paper all high-resolution 
figures are fully zoomable. 

With this initial condition, the tyger phenomenon is par- 
ticularly simple to observe. Fig. Q] shows the solution of the 
truncated Burgers equation with K G = 700 at t = i* and at 
slightly later times Near x = 3tt/2 rj 4.712 the cubic- 
root preshock singularity at t+ and the shock beyond that time 
are standard features. The new one is the "tyger", a growing 
bulge near the positive-strain stagnation point ll38ll . A more 
detailed view is shown in Fig.[2]which shows the discrepancy 



u = v — u between the truncated solution and the untruncated 
one, and zooms into the tyger region and the preshock regions 
at time t+ and shortly thereafter. The bulges seen consist ba- 
sically of oscillations at the Galerkin wavelength 
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with a localized envelope function which is symmetrical 
around the center; this symmetry will actually be destroyed 
later by nonlinearity, as we shall see in Sec. Ill Bl 
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FIG 2: (Color online) Close-up views of the discrepancy u(x) at 
time t = t* = 1 and t = 1.03 in the tyger regions (panels (a) 
and (c), respectively) and in the preshock region (panels (b) and (d), 
respectively). Same conditions as Fig [TJ Note that the tyger bulge 
grows much faster than the preshock bulge. 



We now use a multi-mode initial condition 



uq(x) = sin(x) + sin(2a; + 0.9) + sin(3a;), 



(8) 



for which the first singularity is at t* = 0.2218. Simulat- 
ing the truncated solution, again with K G = 700, we see (cf. 
Fig-E) that at t = 0.25 there is a well-developed shock (near 
x = 1.30), decorated by tygers on each side. 

The tygers are centered around points where the velocity 
equals the half-sum of the limiting velocities when approach- 
ing the shock from the left and the right. For the Burgers equa- 
tion this is precisely the velocity of the shock. Note that there 
is yet another point (around x = 3.30) which has the same 
velocity but no tyger; it has however a negative strain. Thus 
tygers appear to be born at points of positive strain having the 
same velocity as a shock. 

We turn now to the phenomenological explanation of 
the tyger phenomenon, leaving more systematic theory for 
Sec.|nl] The presence of the Galerkin truncation projector P Kg 
makes (0 nonlocal in physical (x) space. A localized strong 
nonlinearity, such as is present at a preshock or a shock, acts 
as a source of a truncation wave whose spatial dependence 
is the Fourier transform of the low -pass filter projector l39ll . 
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FIG. 3: (Color online) Three-mode initial condition. Other parame- 
ters as in Fig.Q] Tygers appear at the points having the same velocity 
as the shock and positive strain. 



More precisely, in physical space the nonlinear term involves 
a convolution with 
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Away from the source (where it is close to a Dirac measure for 
large K G ), this "truncation wave" is mostly a plane wave with 
wavenumber close to K G (see Fig. |4|i and thus has a wave- 
length close to the Galerkin wavelength, 
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Observe that if a preshock/shock is moving with velocity v s , 
the associated truncation wave becomes a progressive wave 
with phase velocity v s . 

Away from the shock region, the Burgers equation can be 
interpreted as describing a particle dynamics: from a La- 
grangian point of view, fluid particles just move with their 
velocity unchanged. In the presence of truncation, those par- 
ticles which happen to have a velocity equal to the phase ve- 
locity of a truncation wave can resonantly interact with such 
waves l40tl . 

Resonant particle-wave interaction is a well known phe- 
nomenon. It is used, for example, in plasma physics to ex- 
plain Landau damping ITill : the near coincidence of the ve- 
locity of particles and of the phase velocity of Langmuir wave 
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FIG. 4: (Color online) Truncation wave with wavenumber K G 
700. 



allows efficient interactions between the two. This can lead 
to wave attenuation (actual Landau damping) or enhancement 
(beam instability, bump on the tail instability, . . . ). There are 
however substantial differences between the Vlasov equation 
(governing the Landau instability) and the truncated Burgers 
equation. For example the Langmuir wave evolves through 
energy transfers via the resonance, a problem which in the lin- 
ear approximation can be solved by use of the Laplace trans- 
form (leading to Landau's rule for a pole-avoiding integration 
path); in contrast, our truncation waves are completely pre- 
scribed by the singularities (preshocks or shocks) and undergo 
no damping in the linear approximation. Furthermore the lin- 
ear approximaton does not have an easy analytic solution (cf. 
Sec. Hill ). In Landau damping, resonant particles get trapped 
and a characteristic "cat eyes" phase-space distribution with 
progressively thinning filamentary structure is obtained (cf. 
Fig. II- 1 of IU9I1 ). In the truncated Burgers dynamics thinning 
is arrested by truncation. 

The radiation of truncation waves begins only at or close 
to the time of formation of a preshock. After a time r has 
elapsed, those fluid particles having a velocity v that does not 
match the preshock velocity v s may not feel much pull from 
the truncation waves if phase mixing is present. More pre- 
cisely, resonant interactions are confined to particles such that 



tAv = t\v 



(11) 



If r is small, as in Figs. [T] and [3j the region of resonance will 
be confined to a small neighborhood of the point of resonance. 
In Sec. lIIIBl we shall see that around the time of appearance 
of a preshock the width of the associated resonance regions 
is typically the order of K G 1 . Outside of such regions, 
including near the preshock, the effect of truncation waves is 
just a small-amplitude oscillation at the Galerkin wavelength, 
which on Figs. Q] and [3] shows up as a thickening of the line 
with respect to the inviscid-limit solution. 

Actually, only resonance points with positive strain produce 
tygers. In a region of negative strain a wave of wavenumber 
close to K G will be squeezed, potentially acquiring a larger 
wavenumber, and thus disappearing beyond the truncation 
horizon which acts as a kind of black hole. 

Observe that in the immediate neighborhood of a preshock 
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or of a nascent shock the strain is also negative and actu- 
ally very large. Although the strongest truncation waves are 
generated near such points (as one would infer by a Gibbs- 
phenomenon argument), their growth is severely hampered by 
the negative strain. Hence the bulge near the preshock grows 
in time much more slowly than that at the positive-strain res- 
onance point, as illustrated in Fig. [2] Actually, the effect of 
truncation near a shock will remain very small (and almost in- 
visible without zooming) until the tyger has fully spread out 
on the ramp (cf. end of Sec. Ill Bb . This situation is in contrast 
with what one observes with other energy-conserving semi- 
discrete schemes |4j[l . such as the dispersive one studied by 
Goodman and Lax 0] . 



B. From tygers to thermalization: the temporal evolution 

It is well known that the Galerkin-truncated solution of the 
Burgers equation, will eventually thermalize to a Gaussian 
state. The simplest case is when the initial velocity vq inte- 
grates to zero over the spatial period. An ergodicity argument, 
supported by numerical simulations, suggests that the ther- 
malized state has equipartition of energy between all Fourier 
modes and thus is just low-pass filtered white noise in the x 
variable id IH. 

How does the highly organized and localized tyger structure 
seen in Figs.[T]and[3]evolve into such a totally random state? 

The birth of tygers around the time of the first preshock 
will be studied in Sec. [HI] In particular, in Sec. lIII Al we shall 
present evidence for scaling properties with K G of both tyger 
amplitude and width at t = t*. Here we shall focus on the 
temporal evolution at later times. Figs. [5] corresponding to the 
single-mode initial condition and K G = 700, show the evolu- 
tion of the tyger in terms of the discrepancy u = v — u from 
t = 1.07 tot = 1.50, shortly after its birth around t* = 1. The 
first few panels display very symmetrical (even) bulges whose 
amplitudes grow in time, because truncation wave input has 
accumulated, while their width decreases (thinning), as a con- 
sequence of phase mixing. Rather quickly, the decrease in 
width leads to a collapse of the tyger, around t = 1.19. This 
is preceded and accompanied by a growing asymmetry of the 
tyger for which we offer the following interpretation. The 
tyger contains kinetic energy in the form of modulated os- 
cillations at the Galerkin wavelength (we shall see that this 
increasing energy compensates the loss of energy in shocks). 
Over scales large compared to the Galerkin wavelength but 
small compared to the tyger width, this kinetic energy gives 
rise to an x-dependent Reynolds stress, which pulls the tyger 
envelope up where the envelope has a negative slope and down 
where the slope is positive. The resulting asymmetry becomes 
very conspicuous after the collapse, as seen in the last panel 
of Fig. [5] This panel has at least two other noteworthy fea- 
tures. To the right and left of the central point x = it we see 
two pieces that look like a portion of (if G -truncated) white 
noise; this is the very beginning of thermalization. Observe 
that the right piece is shifted vertically with respect to the left 
one and that the transition looks almost like an antishock (a 
shock which goes up rather than down, as prescribed by the 



inviscid limit). 

Let us also observe that around the time of collapse there 
is an apparent change of symmetry. Since the single-mode 
initial condition is odd (after shifting the origin of the x-axis 
to the center of the tyger) it stays odd at all later times. This 
is however a statement about the full solution, down to the 
Galerkin wavelength. If we concentrate on the larger-scale 
aspect (the envelope of the bulge), we find that the discrepancy 
is even until somewhat before collapse and odd after collapse. 

For later phases of the tyger growth, it is better to show si- 
multaneously the truncated solution and the untruncated one. 
Also, we switch to the three-mode initial condition which has 
less symmetry and is thus more generic. Fig. [6] shows the 
evolution from t = 0.3, slightly after the first singularity at 
t* = 0.2218, to t = 4.5 when the solution of the truncated 
problem is basically completely thermalized. From t = 0.3 
to t = 0.8 we observe that the chaotic-looking thermalized 
regions born after collapse are growing in extent and are af- 
fecting more and more of the ramp-like structure which are a 
well-known feature of the solution of the untruncated Burgers 
equation after the formation of shocks. In short, we shall say 
that "the tyger spreads out on the ramp". 

As long as significant tyger activity has not reached the 
shocks, their positions, amplitudes and motions are correctly 
described by the Burgers equation, down to the Galerkin 
wavelength. We have checked that during this phase even 
shock merger is unaffected by truncation (an instance is seen 
around t = 1.0). Later, strong tyger activity near the edges of 
shocks is able to shift them slightly (this is visible at t = 1.3). 
Once the shock amplitude has decayed to values much less 
than the tyger fluctuations, the solution looks globally ther- 
malized (t = 4.5). It must be noted that the mechanism 
which prevents thermalization in the Fermi-Pasta-Ulam prob- 
lem 12111 does not seem to be present here. 



C. 2D Euler 

So far we have worked with a very special hydrodynamical 
equation: the Burgers equation is integrable, compressible and 
its solutions generically blow up after a finite time. Is the 
tyger phenomenon also present when none of these properties 
hold, as is the case for the two-dimensional incompressible 
Euler equation with smooth (analytic) initial data and space- 
periodic boundary conditions? 

The short answer is "yes". We have numerically investi- 
gated quite a number of different initial conditions, including 
the two-mode Standard Orthogonal Case (SOC) initial 
condition used in l22ll and random initial conditions. The 
simulations were done with resolutions from 512 2 to 8192 2 . 
Here we report the results for the random initial condition, 
which is 27r-periodic in x\ and X2- The Fourier space 
consists of couples of signed integers k = (ki, fca)- It is here 
decomposed for convenience into shells corresponding to a 
K < |fe| < K+ 1, where K is an integer. Each such shell has 
N(K) Fourier modes. For fc in the K th shell, the Fourier co- 
efficients of the initial vorticity are taken all with the same 
modulus 2K 7 / 2 exp(— K 2 /A)/N{K) and with phases that 
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FIG. 5: (Color online) Evolution of the tyger (discrepancy) for same conditions as in Fig.Q~| growth, thinning, asymmetrization, collapse and 
chaotization. 
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FIG. 6: (Color online) Long-time evolution of three-mode initial condition (other parameters as in Fig.[B. The cyan curve (light grey) is the 
Galerkin-truncated solution and the red one (black) the untruncated inviscid-limit solution. Observe that tygers progressively invade the ramps 
between shocks but that shocks remain sharp and correctly placed as long as the spreading out of the tygers on the ramps has not reached them. 
At very long times, the truncated solution is thermalized. 




FIG. 7: (Color online) (a) shows contours of vorticity with random initial conditions; positive values in dark blue (black), negative values in 
light blue (light gray); (c): same at time t — 0.66 (before appearance of tygers); (b) and (d) are contours of the Laplacian of the vorticity at 
times t = and t = 0.66. Notice the thin elongated red (black) "cigars" which play here the role of the preshocks/shocks in the Burgers case. 
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are uniformly and independently distributed in the interval 
[0, 2tt[, except that opposite wavevectors are given opposite 
phases to preserve Hermitian symmetry. The tyger calcula- 
tions shown here are all with resolution 1024 2 and Galerkin 
truncation wavenumber K G = 342 = (1024 + 2)/3. 
The particular realization used as initial condition 
in the present calculation can be retrieved from 
|http : //www, kyoryu . scphys . kyoto-u . ac . jp/|ht. 
Fig. [7] shows the vorticity and its Laplacian at t = and 
t = 0.66, the latest time at which no tyger is seen (at least in 
the fields displayed). Although for the untruncated solution 
real singularities are ruled out at any finite time, there is 
strong enhancement of spatial derivatives of the vorticity 
ll23ll . The highest values of the Laplacian is found in the 
straight cigar-like structure seen near the center of the figure 
which — as we shall see — will play an important role in tyger 
generation ll42ll . Furthermore this cigar moves very little 
because there is a velocity stagnation point near its center 
HH]- Fig. [8] which is centered on the strongest cigar, shows 
the development of tygers. In terms of the Laplacian of 
the vorticity they become visible around t = 0.71 and then 
become much stronger. 

A further look at a tyger is provided in Fig.|9]which zooms 
into one of the tygers and also shows the Laplacian of the 
vorticity along a cut. 

We immediately see that, as for the ID Burgers equation, 
most of these tygers have come "out of the blue," namely ap- 
pearing at places which had no preexisting small-scale activ- 
ity; more precisely, they appear when complex-space singu- 
larities come within one Galerkin wavelength (2ir) / K G of the 
real domain. The streamlines shown in Fig. [8] indicate that 
tyger activity appears at places where the velocity is roughly 
parallel to the central cigar. As already pointed out, the cigar 
hardly moves; this condition is thus equivalent to having fluid 
particles whose distance to the cigar remains roughly constant. 
In so far as the cigar may be considered as a one-dimensional 
straight object, the truncation waves generated by the cigar 
will have crests parallel to the cigar and those fluid parti- 
cles which move parallel to the crest keep a constant phase 
and thus have resonant interactions with the truncation waves. 
So far, this is basically the same mechanism as discussed in 
Sec. Ill Al for ID Burgers dynamics except, of course, that the 
flow being now incompressible, the velocity within the tyger 
is mostly perpendicular to the direction of fastest variation. 
If we now consider the one-parameter family of straight lines 
perpendicular to a given cigar, each such line will have some 
number (possibly zero) of resonance points; altogether they 
form the tygers. Since the flow outside cigars is fully two- 
dimensional, these tygers have no reason to be parallel to the 
cigars. 

Observe that there are some points where this kind of res- 
onance condition holds but no tyger is seen, for example in 
Fig. [8] at t = 0.71 near x\ = 3.8 and = 3.4. This can 
be interpreted in terms of strain: an incompressible flow has 
at each point a strain matrix with two perpendicular eigendi- 
rections, one for positive strain and the other one for negative 
strain. Fig. [8] has little pink (light gray) segments indicating 
the direction of positive strain. Tyger activity is found only 



at resonance points where the (positive) strain direction is not 
far from being perpendicular to the cigar. More precisely, it 
is easily shown that it has to be within less than ir/4 of this 
direction. Otherwise the near-truncation activity generated by 
resonance is sheared quickly beyond the truncation horizon. 

It is also of interest to show truncation effects and tygers in 
Fourier space. Fig. [TOlshows, at various times, the moduli of 
ketoF^pripxEbeffteifents in the (fci, &2) plane on a logarithmic 
scale. The lowest value contours are at the 10 -15 level, while 
rounding errors are about 10~ 16 . Note that the Fourier-space 
picture is organized in the form of one main lobe, perpendic- 
ular to the physical-space central cigar and secondary lobes 
associated to less intensive small-scale structures. At the ear- 
liest time t = 0.4 no truncation effect near |fc| = K G = 342 
is visible lliill . At t — 0.49, long before tygers become visible 
on the Laplacian of the vorticity in physical space, Fourier- 
space truncation becomes visible. This truncation at first af- 
fects the wavevectors in the direction of the main lobe that 
is perpendicular to the central cigar. Truncation effects then 
spread progressively to other angular directions but appear to 
do so continuously, in contrast to physical space where tygers 
are born "out of the blue." We can actually see such early trun- 
cation effects in physical space by taking more spatial deriva- 
tives and thus putting more weight on high wavenumbers. At 
t = 0.49, Fig. [TT1 shows contours of the tri-Laplacian (V 2 ) 3 w 
with wiggly tygers. By performing various cuts (not shown), 
we checked that the spatial variation is mostly perpendicular 
to the cigar. 

Recently we checked that many features observed for 2D 
incompressible Euler tygers are also present in the 3D case. 
The details will be reported elsewhere. 



D. The dissipative anomaly and the lack of weak limit 

Earlier in this paper, we have seen that for large values of 
the truncation wavenumber K G , the Galerkin-truncated solu- 
tion remarkably preserves many features of the inviscid limit 
such as shocks and their dynamics. So we ask: could it be that 
the Galerkin-truncated Burgers equation converges in a suit- 
able sense to the inviscid limit solution as K G — s- 00? This 
question was actually the main motivation of the present work. 
We shall see that the answer is "no", but a qualified no. 

First, for the kind of analytic initial conditions considered 
here that go singular at some finite time t±, the answer to the 
above questions is actually "yes" for times < t < i*. At 
such times, the solution of the Burgers equation in Fourier 
space for large wavenumbers is bounded by Ce~ Sk (see Ap- 
pendix|B). The effect of truncation is thus exponentially small 
in K G and should go away when K G — > 00. For the two- 
dimensional incompressible Euler equation with periodic an- 
alytic initial data, analyticity holds for arbitrary large times 
and thus the Galerkin-truncated solution is expected to con- 
verge to the untruncated solution l45ll . For three-dimensional 
flow the situation may depend on whether or not there is finite- 
time blow up, a question which is very much open (cf. flolfTlll 
and references therein.) 

Returning to the one-dimensional Burgers case, what about 
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FIG. 8: (Color online) A 2D tyger: before (t = 0.66), early (t — 0.71) and later (t — 0.75). Figures, moderately zoomed, centered on the 
main cigar. Contours of the Laplacian of vorticity in red (black), ranging from —200 to 200 by increments of 25, streamlines in gray, ranging 
from —1.6 to 1.6 by increments of 2 and positive strain eigendirections in pink (light gray) segments. Notice tygers appearing at t — 0.71 
in the form of wavy red (black) contours at places where the velocity is roughly parallel to the main cigar; later, the tygers grow in strength 
and extension. A blue (black) x-mark is added at t = 0.71 to underline that there is no tyger, in spite of resonance, because of a wrong strain 
direction. 
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FIG. 9: (Color online) Upper: zoomed version of contours of the 
Laplacian of vorticity at t = 0.71. Lower: plot of the Laplacian of 
vorticity along the horizontal segment near X2 = 3, shown in the 
upper panel. 



t > t*, when shocks are present and the solution is dissipa- 
tive, whereas the Galerkin-truncated solution is conservative! 
How can such a conservative system mimic the dissipative 
anomaly? One could imagine that the small-scale tyger ac- 
tivity plays the role of molecular motion and that motion on 
scales much larger than the Galerkin wavelength is governed 
by the inviscid-limit Burgers equation. 

This is however not the case. Tadmor investigated the 
limit — in a suitable sense to be defined below — of the 
truncated solution and found that it cannot be dissipative 
(Ref. iflill p. 31). The limit considered by Tadmor and in 
other papers studying conservative modifications of the Burg- 
ers equation with strong oscillations is a distributional 
weak limit. Assuming that the Galerkin-truncated solution v 
has a weak limit v satisfying the Burgers equation, and us- 
ing the basic dynamical equations ©-(O, Tadmor shows that 



has the weak limit 



from which he infers that the 



limit is actually a strong one which implies energy conserva- 
tion and contradicts the dissipative character of the solution to 
the Burgers equation. Recently numerical simulations of the 
inviscid Galerkin-truncated Burgers equation with K G up to 
about 10 4 showed indeed that such solutions do not converge 
to the inviscid limit of the untruncated solution l25ll . 

The simplest instance of such a weak limit is just to apply 
a low-pass filter to the solution with a fixed threshold K for 
the modulus of the wavenumber, while letting K G —> oo. It is 



FIG. 10: (Color online) Contours of the modulus of the vor- 
ticity Fourier coefficients at various times. Negative ki values 
not shown because of Hermitian symmetry. Contour values are 
10 _1 , 10 -2 , . . . , 10 -15 from inner to outer (in color version, green, 
blue and pink highlight the values 10 -5 , 10 -10 , and 10~ 15 , respec- 
tively). Galerkin truncation effects are visible above the rounding 
level already at t = 0.49 and become more and more invasive. 




FIG. 11: (Color online) Contours of tri-Laplacian of the vorticity 
(upper: zoomed; lower: superzoomed) showing a tyger already at 
t = 0.49. 
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FIG. 12: (Color online) Plots of solution of the Galerkin-truncated 
Burgers equation, with Kg = 5, 461 in cyan (light grey) and 
Kg = 21, 845 in black, low-pass filtered at wavenumber K = 100, 
at various times. Initial condition vo (x) = sin(x) +sin(2a; — 0.741). 
The untruncated solution is shown in red (black, mostly straight line). 



then easy to show that one can find a subsequence 711,712, .. . 
of the sequence of integers K G such that the low-pass filtered 
solution has a limit 12611 . We have obtained evidence that, 
without taking subsequences, there is no weak limit. Indeed 
Fig.[T2lshows the low-passed solution with the same threshold 
K = 100 for the same initial condition and the same output 
time, but for two very large and well-separated values of the 
truncation: K G = 5,461 (denoted 5K) and K G = 21,845 
(denoted 2 IK). The solutions agree very well with the untrun- 
cated solution at shocks and nearby but the 5K and 21K tygers 
differ significantly, even after application of the low-pass fil- 
ter. It may thus be that there is no weak limit as K G — > 00 
& 

We can supplement this by a more physical and fluid me- 
chanical explanation of why the truncated solution cannot 
converge (weakly) to the inviscid-limit solution. We have seen 
in Sec. lIIBI that shocks behave just as predicted by the invis- 
cid limit for a substantial length of time (until tyger spreading 
on the ramp reaches the shocks). Because of this, the hypo- 
thetical limit of the truncated solution would be losing energy 
at the shocks, just as the ordinary Burgers equation. The en- 
ergy lost has to be found in the tygers in the form of high- 
wavenumber oscillations. Decomposing the truncated solu- 
tion v — u + u (where u is t he inv iscid-limit solution) we ob- 
tain tyger Reynolds stresses u 2 /2, where the overline means a 
mesoscopic spatial averaging over a distance large compared 
to the Galerkin wavelength and small compared to any macro- 
scopic scale (for example using a low pass filter with threshold 
K = 100 for the case K G = 5, 461). If the mesoscopic tyger 
energy and thus the Reynolds stress is not spatially uniform, 
its gradient will drive the flow away from the inviscid limit. 
This is the same mechanism that makes the tyger asymmetri- 
cal, as already mentioned. Is there a way to obtain the correct 
inviscid limit by eliminating the undesirable Reynolds stresses 
through some kind of tyger purging? We shall come back to 
this important practical issue in the last section. 

Finally, let us remark on our choice of the word tyger for 
the oscillations which are a result of Galerkin truncation. His- 
torically, the distinction between conservative and dissipative 
systems has played a crucial role in not only man's scientific 
pursuits, but also in a deeper cultural context. For centuries 
there was a certain sanctity associated with things conserva- 
tive as opposed to being dissipative. Hence, before Galileo's 
telescope revealed the "transient" nature of celestial occur- 
rences (e.g., Sun spots), man had always ascribed heavenly 
objects as conservative and the more mundane, transient, and 
earthly occurrences as dissipative. The fine balance of the two 
seemed essential for all existence. 

In the backdrop of this and in our investi- 
gations of truncated systems, William Blake's 
poem "The Tyger" assumes special significance 
|http : //en . wikipedia .org/ wiki/The_Tyger| 
In most interpretations, Blake's tyger is of course not the 
animal itself |H3l; it is a metaphor for the symmetry of 
seemingly different, even opposite, processes which nonethe- 
less combine to make a coherent whole. Thus apparent 
oppositions such as life and death, light and darkness are 
seamlessly unified. In our present study, we explore the 
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interplay of conservative and dissipative dynamics: how one, 
surprisingly, might be embedded into another. As a result, it 
was quite natural for us to call this phenomenon a tyger. 

IIL DETAILED ANALYSIS OF THE BIRTH OF THE 
BURGERS TYGER 

So far our point of view has been that of the (numerical) ex- 
perimentalist, with some amount of phenomenological theory 
used to interpret the results whenever possible. In this section 
we shall use a lot more "systematic theory", namely expan- 
sions for large values of the Galerkin truncation wavenumber 
K G . In principle, such expansions can be carried out beyond 
leading order; we shall however not attempt this, let alone ob- 
tain rigorous bounds for errors. Indeed some of our approxi- 
mations used below are akin to what one calls "patching" (as 
opposed to "matching") in boundary layer theory, that is ap- 
proximations which can give the correct exponent of a power- 
law leading term but cannot predict the correct constant in 
front. 

Let us now give a general overview of how we intend to 
proceed. We shall concentrate on a single problem, that of 
the birth of the tyger for the one-dimensional Burgers equa- 
tion with single-mode initial condition. The birth, which takes 
place around the "Galerkin time" t G when complex singular- 
ities come within one Galerkin wavelength from the real do- 
main, is part of an early phase which extends from t = to 
the time of the first singularity t = t* = 1. In Sec. IIII Al we 
shall see that the early tyger at t = has remarkable scaling 
properties with K G . Our intention here is to understand ana- 
lytically how this comes about. For this, our strategy will be to 
devise various models/approximations which make the prob- 
lem simpler while hopefully keeping the leading-order behav- 
ior unaffected. There will be three levels of modelling: (i) 
linearization around the untruncated solution, (ii) ignoring the 
"exponentially small" phase up to the Galerkin time t G , and 
(iii) "freezing", i.e. replacing the untruncated solution u(x, t) 
by u+(x) = u(x, t±) for t e [t G , t*]. In Sec. IIII Bl we explain 
how to do this and why it is justified. The problem is then re- 
duced to studying a linear first-order differential equation with 
constant coefficients in a finite dimensional space (Sec. IIII C| >. 
It was brought to our attention recently by J. Goodman that 
this differential equation has features in common with that 
studied by Goodman, Hou and Tadmor (subsequently cited 
as GHT) in connection with the stability of the pseudospectral 
method in the presence of aliasing 12811 . 



A. Scaling properties of the early tyger 

In Sec.HJwe nave seen that tygers are born at suitable res- 
onance points in the form of bulges made of oscillation at the 
Galerkin wavelength with a very symmetrical envelope. Even- 
tually, Reynolds stresses will distort this envelope. At the time 
t± of the first singularity this is not yet the case. Now we con- 
centrate on the scaling properties for high K G of the amplitude 
a and width w of such early tygers, in terms of the discrepancy 



u = v — u. Here we limit ourselves to the single-mode initial 
condition. Figs. [T3l and [T4l show respectively the amplitude 
and width as a function of K G for values ranging from 100 to 
40,000. The amplitude is measured at the maximum closest 
to the center of the tyger, located a distance it/ (2K g ) away. 
The width is measured at half of this amplitude. It is seen that 
both have clean scaling laws for sufficiently large K G : sub- 
dominant corrections are perceptible on the log-log plots only 
fori^c < 200. The leading-order behavior is 

acxK G - 2/3 , w(xK G - 1/3 . (12) 

The following sections are mostly devoted to explaining these 
scaling laws. 
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FIG. 13: (Color online) A log-log plot of the discrepancy-based am- 
plitude a of the tyger at time t — t* as a function of A'g. Initial 
condition uo (x) = sin x. The data points are shown by filled red 
circles (black circles) and the thick black line is the best power-law 
fit oc Kq 2 ^ 3 , which holds over two decades. 
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FIG. 14: (Color online) Log-log plot of the width (at half amplitude) 
of the tyger. Otherwise, same as in Fig.[T3] The best power-law fit is 
now oc Kq 1 ^ 3 , which holds over nearly two decades. 



B. Three successive approximations 

As already observed in Sec. Ill Dl initial conditions with a 
finite number of Fourier harmonics, such as the single-mode 
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case, are for a while analytic in the complexified space vari- 
able within a strip around the real domain of width 6(t). As 
long as 6(t)K G 1 the effect of truncation is exponentially 
small. This is why, the kind of tygers reported in Sec. [II] are 
not seen before t G , when complex singularities come within 
roughly one Galerkin wavelength from the real domain. For 
large K G , this happens only a short time 0(K G ~ 2 / 3 ) before 
the time t* of the first singularity (cf. Appendix [B). Hence, 
by the time t*, truncation has been felt significantly only for a 
lapse of time 0(K G ~ 2 / 3 ). The phase mixing argument given 
in Sec. Ill Al and in particular (fTTT i tell us that the coherent build 
up of a tyger will affect only those locations whose velocity 
differs from that at resonance by an amount Av such that 

^i^^" 3 <13) 

Since at such times, the velocity v of the truncated solution 
is expected to stay close to the velocity u of the untruncated 
solution and the latter varies linearly with x near the reso- 
nance point, the width of the tyger is itself proportional to 
K G 1 , as indicated by the simulations. 

Is there an equally simple argument to understand why the 
amplitude scales as A' G ~ 2 / 3 ? One possibility would be to ob- 
serve that at the untruncated solution has a cubic root be- 
havior near the preshock location (cf. Appendix IBl. If we cut 
out a small interval of length one Galerkin wavelength A G , 
we will be "missing" an energy ~ J G x 2 l 3 dx ~ K G ~ 5 ^ 3 . 
Remembering that the Galerkin-truncated equation conserves 
energy, if we assume that this missing energy is transferred 
entirely to a tyger of width A'q" 1 / 3 , we obtain precisely an 
amplitude oc K G ~ 2 ^ 3 . There is however no reason to assume 
that near the preshock the effect of truncation can be reduced 
to carving out a little interval of one Galerkin wavelength. It 
is doubtful that this energy argument can be turned into some- 
thing rigorous. 

We turn now to more systematic arguments. Before i*, 
when S(t) > 0, it is easily shown that the Galerkin-truncated 
solution converges strongly to the untruncated solution. The 
simulations reported in Sec. IIIIAl suggest that this still holds 
at t = i*, since the amplitude of the discrepancy goes to zero. 
Furthermore, the nonlinear effect of Reynolds stresses, as dis- 
cussed in Sec. Ill Dl is of even higher order. Indeed, within 
the bulge the Reynolds stresses will be ~ K G ~ 4 '' 3 ; since 
they change spatially on a scale ~ K G ' , the gradient of 
Reynolds stresses is ~ A" G _1 ; over a time interval ~ K G ~ 2 ^ 3 
this will change the bulge amplitude by ~ K G ~ 5 ^ 3 , which is 
small compared to the amplitude of the bulge itself. 

All this suggests that the early tyger development can be 
captured by somehow linearizing the Galerkin truncated so- 
lution around the untruncated one. Let us rewrite the basic 
dynamical equations ©-(O in terms of the discrepancy 



u = v — u. 



(14) 



We obtain 

dtu + P KQ d x (uu + ±\=(l- P Kg ) x IL fi(o) = 0, 

(15) 



where I stands for the identity operator and the zero initial 
condition follows from uo = Vq, a consequence of having a 
finite number of modes initially. 

Observe that the r.h.s. of (fT~5T > provides no input to 
wavenumber below the truncation. Actually this input is hid- 
den in the l.h.s. To make this clear, we need to decompose 
the various fields into their Galerkin-truncated part and the re- 
mainder. We set 



E. u 

K G 



U> EE (I - R W 

V K G ) 



(16) 
(17) 



Next, to similarly decompose u, we use (I — P 1<g )u 



which follows from the fact that v = u 
beyond the truncation. Hence we have 



u = u — u 

u' EE E, U. 

K G 



u has no harmonics 

(18) 
(19) 



In what follows we shall work mostly with v! which has no 
harmonics beyond the truncation and which we call the per- 
turbation. As we shall see, around t*, the perturbation is 
small. In contrast, beyond the truncation, the discrepancy is 
just equal to minus the untruncated flow; it is thus known but 
in no way small. Now, we apply P^ g to dT3T > and use the vari- 
ous decompositions to obtain 



d t u'+P K d x 



W) 2 



P K d x u<u> 



(u>f 



(20) 

The r.h.s. of (120b is a known function which we shall call the 
beating input and denote / because it describes how harmon- 
ics of the untruncated solution, located beyond the truncation, 
interact with themselves or with subtruncation harmonics to 
give a subtruncation input. This beating input, which is shown 
in the lower panel of Fig.[l9j consists basically of spatial os- 
cillations at the Galerkin wavelength, modulated by an enve- 
lope that peaks at the preshock. This is the precise content of 
what we called "truncation waves" in the phenomenological 
approach of Sec. Ill Al 

We are now in a position to define the three approximations 
made for large K G , which we call respectively linearization, 
reinitialisation and freezing: 

• 1. The term (u') 2 in the l.h.s. of ( f20T > is discarded; 

• 2. The perturbation vl is set to zero at time t G \ 

• 3. The untruncated solution is frozen to its value. 

Concerning linearization, we already observed that nonlin- 
ear effects will be weak if the scaling laws for the i* tyger 
are indeed given by (fT2l . Reinitialization is justified because, 
prior to the Galerkin time t G , only exponentially small per- 
turbation are present and, here, we are only trying to capture 
algebraically small terms. It is then convenient to introduce a 
new shifted time variable 



t - tr 



(21) 
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To avoid unnecessary constants, we choose, t G = t* — 
K G ~ 2 ^ 3 . Hence, in the t variable, the first real (preshock) 
singularity is at 



simply given by 



(22) 



We shall also take the liberty to still denote the perturbation 
by v! when it is expressed in terms of the shifted time. As to 
the freezing, replacing u(t) by = u(t+), it appears justi- 
fied since the untruncated solution hardly changes between t G 
and t+, except at wavenumbers much larger than K G which 
do not contribute to (1201 1. Freezing changes (|2Q| i into an equa- 
tion with constant coefficients. In particular the beating input 
becomes time-independent. 

Although there are good theoretical reasons to make 
these three approximations, we also tested them numerically. 
Fig. [15] shows the scaling law for the tyger amplitude at t+ 
in compensated form (after multiplication by K G +2 ^ 3 ) for (i) 
the full problem, (ii) with only linearization assumed and (iii) 
with in addition reinitialization and freezing assumed. All 
three cases have the same scaling law. Linearization brings 
about only a minuscule change, as expected. The other two 
approximations make a difference of about twenty percent, an 
indication that they can be refined. 
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FIG. 15: (Color online) Log-linear plot of the compensated ampli- 
tude of the tyger Kq^ 3 a(K G ), calculated (i) from ft, (ii) from the 
linearised approximation for u', and (iii) from the freezing plus reini- 
tialization approximation, all versus K G . The data corresponding 
to (i) and (ii) are indistinguishable and shown as blue filled circles 
(black upper curve). The data corresponding to (iii), shown by black 
filled circles with a red border (black lower curve), have the same 
scaling but a multiplicative constant about twenty per cent lower. 



C. Tyger birth: the minimal model 

With all three approximations formulated in the previous 
section, the temporal dynamics of the perturbation near is 



Au + /, it'(0) = 0, 



dr 

A = -P K dx («*•) , 



f = p KG d x [utu> 



«) 2 



(23) 
(24) 

(25) 



It is important to stress that, because of the freezing ap- 
proximation, A and / are evaluated at t* and thus time- 
independent. 

We shall now write this equation as a finite-dimensional lin- 
ear differential equation by working in Fourier space and see 
that A is actually a matrix. Some of the quantities that we shall 
look at are more conveniently represented by taking the initial 
condition uq = sin x, which has the tyger born near x = 
and the preshock at x = it. For other quantities the choice 
uq = — sin x is better. We shall refer to the former as "origin 
at the tyger" and to the latter as "origin at the preshock". Note 
that in both cases, the initial condition being real and odd, so 
are the untruncated solution u, the truncated one v, the pertur- 
bation u' and the beating input /. Hence their Fourier coeffi- 
cients Uk, Vk, u' k and //. are pure imaginary and odd functions 
of k. With this notation, it is easy to rewrite d23l-(l25T) as a 
system of 2K G + 1 equations, indexed by k E 



-K G , K G }: 



Ac: 



-t-K = Akk ' "'fe' + h 

k' = -K G 



A kk > 



-ik u± 



k-k> 



f k = ik ^ i^p u+q + -jKp 



i' fc (0) = 0,(26) 

(27) 
(28) 



p+q=k 



It is important to observe that, u< being pure imaginary, 
the entries of the matrix A are all real. 

Here we observe that GHT were led, at the technical level, 
to studying a homogeneous version of d26| i without the beat- 
ing input / (and thus without resonant wave interactions). In 
their work, the velocity u is also prescribed but mostly taken 
to be sin(x) or s'm(px). The operator/matrix Akk' differs only 
marginally from ours, due to the deliberate presence of alias- 
ing. 

It is of course quite easy to solve j26l)-(l28l) numerically for 
the perturbation u'(t*). Fig. [TBI shows its imaginary part for 
three large values of K G . Most of the activity is concentrated 
in boundary layers near ±K G (we are only showing the right 
boundary layer because the function is odd). In the "tail" out- 
side of this boundary layer the amplitude of the perturbation 
is very small. The Fourier space boundary layer contains all 
the information about spatial oscillations on scales close to 
the Galerkin wavelength. If a tyger is localized in physical 
space near X, the Fourier amplitude will have a phase factor 
e~ lkx . For this particular set of figures we have taken the ori- 
gin at the preshock. In the boundary layer the most conspic- 
uous features are even-odd oscillations which are a signature 
of a tyger located at X = it. The oscillations are not com- 
pletely symmetrical between positive and negative values, an 
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FIG. 16: (Color online) The boundary layer in Fourier space near K G . Shown are the imaginary parts of u'(i*) for three values of K G , The 
origin is at the preshock. The even-odd oscillations indicate that most of the activity is at the tyger, a distance -n away. 



indication that there is also some small-scale activity near the 
preshock at X = 0. The low-amplitude tail to the left of the 
boundary layer has no oscillations, indicating that it comes 
mostly from the neighborhood of the preshock. The width of 
the boundary layer is found to scale approximately as K G ' , 
that is the inverse of the tyger width, as expected from Heisen- 
berg's "uncertainty principle" f48f|. As to the peak amplitude 
of the perturbation in the boundary layer, it is found to scale 
as K G ~ . After subtraction of the high-fc contribution stem- 
ming from the neighborhood of the preshock j49tl . we found 
that the upper part of the envelope of the boundary layer has 
the following scaling representation: 

Im^ = K^F (j^TJt) ' F ^ = O.U8e~ 3A5x , 

(29) 

which implies a collapse of all the boundary layer data after 
suitable rescaling, as illustrated in Fig.fTTl 
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FIG. 17: (Color online) The envelopes of the various boundary lay- 
ers shown in Fig. [16] (with preshock contributions subtracted out), 
collapsed into a single curve after rescaling. Red (black) circles: 
K G = 20, 000, blue (black) squares: K G = 15, 000, red (black) 
triangles: K G = 10, 000, blue (black) diamonds: K G = 5, 000. The 
thick black line is the exponential fit (f29). 

Note that by Fourier transformation, (|29l immediately im- 
plies the basic scaling laws (fTSl i for the width and amplitude 



of the tyger. As we have seen, the K G ' dependence of 
the width is just a consequence of the phase mixing condition 
(fTTT i. In the next section we shall try to understand how the 
amplitude factor comes about. 

D. Tyger birth reduced to a simple (?) linear algebra problem 
Obviously, ( f23l > can be solved for the perturbation at time 

7* = ^g" 

«'(t*) = A" 1 (e^ -I) / = f Y^A-^^J f. (30) 

If the matrix A is singular (as it actually is), the middle equa- 
tion of (f30b is not directly meaningful, but the r.h.s. remains 
meaningful. 

From this it becomes clear that much is controlled by the 
spectral properties of the operator A. We are thus led to con- 
sider the associated eigenvalue/eigenvector equation 

Aip = Xip, (31) 

which plays for the Galerkin-truncated problem the role of 
the standard Orr-Sommerfeld equation |221 and will be thus 
called. Detailed spectral properties of the Orr-Sommerfeld 
operator A, which is neither Hermitian nor antihermitian, are 
discussed in Appendix IDl The eigenvalues, most of which 
are complex, come in opposite pairs; the associated complex 
eigenvectors being either even or odd functions of k. In addi- 
tion there is a zero mode, i.e. an eigenvector with eigenvalue 
zero. We denote the eigenvalues by Xj where j is a signed 
integer varying from — K G to K G . For positive j, Xj is the 
jth eigenvalue with positive imaginary part, eigenvalues be- 
ing ordered by increasing moduli; A_j = — Xj and A = 0. A 

complete set of complex eigenmodes is denoted ip^> and their 
(7) 

Fourier coefficients by ipl . When K G is large, the eigenval- 
ues with large indices j are almost pure imaginary and the 
largest (in moduli) eigenvalues are very close to HK G l50ll . 
The typical spacing between the moduli of successive eigen- 
values is order unity for large K G . 
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Also of importance are the strength and scaling properties 
of the beating input, which are discussed in Appendix [C] Its 
Fourier coefficients ff. are pure imaginary and odd functions 
of k. We saw that the tyger has the unexpected feature that 
it appears away from the preshock. However the beating in- 
put, which is caused directly by truncation, is mostly localized 
where the untruncated solution displays its highest small-scale 
activity, namely near the preshock. With the origin taken at 
the preshock, the imaginary part of the beating input Im f/. 
peaks at truncation k = K G with \f Ka \ = O ^A' G ~ 1/3 ^. 

When moving down from K G to lower k, it falls off rather 
slowly as (K G — fc) -1 / 3 while keeping a constant sign. Thus 
in Fourier space we have a rather broadband beating input. In 
the physical space the beating input is mostly an oscillation at 
the Galerkin wavelength, whose envelope falls off as |a;|~ 2 / 3 , 
where \x\ is the distance from the preshock location. 

Let us now decompose the beating input / and the solution 
m'(t*) in terms of the eigenmodes: 

U '( T *) = £ f = J2 /°'¥ j) - (32) 

3 3 

It then follows from ( T30b that 

u /( %*) = — f {j \ (33) 

A 3 

in which it is understood that the fraction takes the value r* 
when A = 0. An important role is played by those eigenvalues 
for which 

|A j -|t*~1; (34) 

these will be called threshold eigenvalues. Since t* oc 
K G ~ 2/3 there is a whole range of eigenvalues well below and 
well above the threshold. The corresponding eigenmodes will 
be called "low-lying" and "high-lying" modes, respectively. 
Well below threshold the fraction in d33l can be Taylor ex- 
panded, yielding to leading order r*/^: those modes are es- 
sentially unaffected by their interaction with the flow through 
the Orr-Sommerfeld operator A. We have checked that such 
modes are responsible for the low-amplitude oscillatory tail 
to the left of the boundary layer seen in Fig. [16] Well above 
threshold, the modulus of the fraction is much smaller than 
because c AjT * is essentially a phase factor lf5lll and one can 
thus suspect that high-lying modes do not contribute much to 
the boundary layer. Making this argument solid requires a bet- 
ter control over the phases than we have been able to achieve 
analytically. We have thus carried out numerically partial 
summations of the r.h.s. of d32l ), on the one hand starting from 
low-lying eigenmodes and adding progressively higher-lying 
eigenmodes and on the other hand doing it in reverse. In both 
instances we found that the boundary layer of Fig.[T6lemerges 
mostly from modes near the threshold 1521 . 

When d33l is applied near threshold, the 1/Aj yields a fac- 
tor t* oc K G ~ 2 ^ 3 , while the beating input has an overall factor 
oc K G ~ X ^. Together, this produces a I\ G 1 amplitude factor 
in a boundary layer of thickness oc Kq 1 ' 3 , needed to explain 



the K G ~ 2 / 3 law for the amplitude of the t — tyger in physi- 
cal space. What we just explained is however far from a proof 
since (i) we did not show analytically that the dominant con- 
tribution to the boundary layer comes from threshold modes 
and (ii) a K G _1 amplitude factor for the beating input / does 
not necessarily imply the same factor for its threshold compo- 
nents /W. 



IV. OPEN PROBLEMS AND CONCLUSION 

We must now conclude our adventures in Tygerland. Al- 
though this project has been unfolding over three years, we 
have the feeling that we only indented the subject, as far as 
true mathematical understanding is concerned. For exam- 
ple in Sec. [TTTJ devoted just to the birth of tygers, we have 
not identified analytically the important mechanism which al- 
lows threshold modes to populate the boundary layer seen in 
Fig.[0>] 

As to the after-birth events, they have so far only been 
the subject of numerical experimentation and occasional phe- 
nomenological theory. The collapse of the tyger, shortly after 
the time of appearance of the first shock, seen in Fig. [5] is 
strongly reminiscent of the collapse phenomenon in plasmas 
ll3lll . The immediately subsequent development, as we have 
seen in Sec. Ill Bl involves at least two phenomena. One is 
the spreading out on the ramp, the moving of the tyger along 
the ramps of the untruncated Burgers solution, which perhaps 
can be explained by advection effects; the other one is that the 
small-scale motion looses the highly organized structure seen 
around t* ; in other words the Fourier spectrum broadens away 
from the Galerkin wavenumber. This may signal the onset of 
the thermalization of the solution which, eventually, becomes 
a Gaussian noise in the space variable with a flat spectrum. 
Here a digression is in order. The ordinary untruncated Burg- 
ers equation — with or without viscosity — has played a ma- 
jor role, not only as testing ground for numerical schemes, 
but also for helping us to find mistakes in excessively naive 
ideas intended for Navier-Stokes turbulence. The Galerkin- 
truncated Burgers equation may take us a step further, being 
paradoxically closer to Euler-Navier-Stokes: it is nonlocal (in 
a way consistent with energy conservation) and its solutions 
display spatio-temporal chaos, as documented, e.g., in lfl5ll . 

Of course all this would be quite academic if we did not al- 
ready have good evidence that the key phenomena associated 
to tygers are also present in the two-dimensional incompress- 
ible Euler equation, as discussed in Sec. Ill CI Our understand- 
ing of this, so far, based on what we know about the analytic 
structure of 2D flow, is far from complete. It seems important 
to find out how the more systematic theory of the birth can be 
carried over to 2D — and perhaps 3D — Euler. So far, it is clear 
that complex-space singularities approaching the real domain 
within one Galerkin wavelength are the triggering factor, as in 
the ID Burgers case. 

Now, a few remarks about tyger purging, which is defi- 
nitely not the central issue of the present investigation. We 
have seen in Sec. [TT] that tygers, being born far from shocks, 
do not modify shock dynamics but do modify the flow else- 
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where because the tygers induce Reynolds stresses on scales 
much larger than the Galerkin wavelength; hence the weak 
limit of the Galerkin-truncated solution as K G — > oo is defi- 
nitely not the inviscid limit of the untruncated solution. Can 
we "purge tygers away" and thereby obtain a subgrid-scale 
method which describes the inviscid-limit solution right down 
to the Galerkin wavelength? 

How to do this practically? Several ideas come to mind. 
One is simply to apply some amount of viscosity. If the 
viscosity is sufficiently large, the truncation becomes ir- 
relevant but the solution thus obtained will not coincide 
with the inviscid limit down to the Galerkin wavelength or 
anywhere close to it. A second idea is to look for embryonic 
tygers in physical space and selectively abort them. This can 
perhaps be done by a suitable wavelet or filtering technique 
but may be tricky 15311 . A simpler idea is to purge the 
boundary layer near K G at each time step. However this 
amounts to applying a Galerkin truncation with a slightly 
smaller K G and will produce more tygers. A more subtle 
way worth exploring is to wait until a low amplitude tyger 
has appeared that is concentrated in a sufficiently narrow 
boundary layer near K G ifsill and then to perform the purging, 
an operation which clearly should not take place too often. 
Such ideas will of course have to be tested carefully in future 
work. One may also wonder to what extent such a purging 
technique can be carried over to 2D and 3D incompress- 
ible flow. We have already checked that for the case of 
2D and 3D incompressible flow, the birth of the tyger takes 
place in a narrow boundary layer near the Galerkin truncation. 
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Appendices 

Appendix A: Numerics and graphical representation 

With the exception of the inviscid-limit solutions of the un- 
truncated Burgers equation which were obtained by the Fast 
Legendre Transform method 13311 . all the numerical simula- 
tions presented in this paper used the pseudo-spectral method 



Jgt], in combination with a fourth-order Runge-Kutta time 
marching with a time step 10~ 4 for K G < 5, 000 and 10 -5 
for K G > 5, 000. Up to 2 16 collocation points were used. In 
order to implement Galerkin truncation it was essential to re- 
move aliasing. In principle this can be done by the so-called 
two-thirds rule: K G < (2/3)if roax = N/3, where N is the 
number of collocation points. This allows a calculation of the 
solution at collocation points which suffers only from double- 
precision rounding errors and temporal truncation. 

However, when it comes to representing tygers graphically, 
using the two-thirds rule can produce a stroboscopic graph- 
ical artefact, illustrated in Fig. [18] Since in the presence of 
tygers there is a strong excitation at and near the truncation 
wavenumber K G , the velocity is very close to a sine wave with 
Galerkin wavelength. Unless proper interpolation is used, 
such a sine wave cannot be correctly represented using only 
three points per wavelength. Otherwise the rapid oscillations 
disappear in favor of an illusory triple valuedness. One easy 
way to do the interpolation is to use a much higher number of 
collocation points whenever graphical output is needed. 
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x 

FIG. 18: (Color online) Plots of the solution v(x) of the Galerkin- 
truncated Burgers equation, with A' G = 85, using two different num- 
bers of collocation points, namely, N — 16, 384 (black continuous 
line) and with red-filled circles (black circles) for N = 256 for 
the single-mode initial condition vo(x) = sin(x — tt/2) shown at 
time t — 1.10. This plot illustrates the stroboscopic effect, a purely 
graphical interpolation artefact, caused by insufficient resolution for 
a given K G . 



Appendix B: Real and complex singularities for the Burgers 
equation 

As is well known, the solutions of hydrodynamical equa- 
tions such as the Burgers or Euler equations in any dimen- 
sion with space-periodic and analytic initial data, remain so 
for at least a finite time (cf., e.g., Refs. iflolfTlll and references 
therein). In the ID Burgers case, a real singularity appears af- 
ter a finite time (finite-time blow up). For 2D Euler, analyticity 
is preserved for all times but complex-space singularities will 
typically approach the real domain arbitrarily close at long 
times. For 3D Euler it is not clear if there is or not finite-time 
blow up. The signature of complex-space singularities — more 
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precisely of those closest to the real domain — is an exponen- 
tial fall-off at high wavenumbers k of the spatial Fourier trans- 
form oc e _(5 w fe (within algebraic prefactors), where S(t) is the 
width of the analyticity strip, that is the distance from the real 
domain of the nearest complex-space singularity. As long as 
6(t)K G 3> 1 the effect of Galerkin truncation is exponentially 
small and can be ignored for all the purposes of the present 
study. For convenience we define the threshold as the time t G 
when S(t G )K G = 1. 

Let us here recall how S(t) can be obtained explicitly for 
the case of the untruncated Burgers equation with the initial 
condition uq(x) = — sin a; (for details, cf. Ref. ifl^ ). In La- 
grangian coordinates, the motion of an inviscid Burgers fluid 
particle is given by 



imaginary, odd and given to leading order at large wavenum- 
bers by 



(C2) 



where C is a real constant and it is understood that fc~ 4 / 3 has 
the same sign as k. Actually, dC2| ) is a very good representa- 
tion of the Fourier transform of u k at all but a few low-lying 
wavenumbers. Observe that the r.h.s. of ( IC1I ) has two terms, 
the first has \p\ < K G and \q\ > K G while the second has 
both \p\ and \q\ greater than K G . The latter is easily bounded 



in modulus by D\k\K G 



-5/3 



where D is a positive constant 



and thus will be seen to contribute negligibly to the asymp- 
totics. Hence, from dCll ) and dC2l ), we have 



a i— > x = a + tuQ (a) = a — t sin a. 



(Bl) 



A singularity is obtained when the Jacobian dx/da of this 
map vanishes. For t slightly less than = 1, the Lagrangian 
and Eulerian locations of the singularity nearest to the real 
domain is close to the origin and can be obtained by expanding 
sin a to cubic order. We thus have, up to higher-order terms: 

3 2 

x = (t* - t)a + — ; dx/da = t* - t + — . (B2) 
6 2 



The Jacobian is seen to vanish at a*(t) = ±i^/2(t* — t). The 
corresponding Eulerian singularities are at 



fk 
9k 



-C 2 \kg k 

E 

p+q=k, \p\<K G , \q\>K G 



(C3) 

p-^q"^. (C4) 



We introduce dimensionless variables: 



g 

K G 



and rewrite dC4t as 



9k = -K a -*/ 3 p- 4/3 (l+P-k) 



-4/3 



fc<p<l 



(C5) 



(C6) 



(B3) 



Equating the modulus S(t) of the imaginary part of x±(t) to 
1/K G , we obtain: 

/ 3 \ 2/3 

U-t G = i^=\ K G - 2 ^ » 1.04X G - 2 / 3 . (B4) 

For convenience we have replaced 1.04 by unity and thus used 
in Sec. |m] 



i — i — K ~ 2 / 3 

<<G — t'* 



(B5) 



Appendix C: The beating input 

Our purpose it to find the large-_ftT G asymptotic behavior of 
the beating input d28| i, repeated here for convenience, 



h 



k = ^ ( 

p+q=k 



1 



(CI) 



in the range 1 <C K G — k <C K G . Although it is this input 
which eventually permits the birth of the tyger, the beating in- 
put is strongest at the preshock and it is best to place the origin 
there. The function u*, whose high- and low-passed filtered 
Fourier transforms appear in (ICll l, is odd and has a cubic root 
singularity at the origin. Hence its Fourier transform is pure 



where the summation on p is over integer multiples of 1 /K G . 
Now we set 



y 



p — k 



(C7) 



and approximate the summation by an integral to obtain, to 
leading order, 



9k 



- l -K G -V\K a -k)-V\ 



(C8) 



whence 



1 



J k ~ i-C*K Q - y i*(K Q - k) 



K K G - k < K G . 

(C9) 

There are also subleading corrections involving (K G — k)°, 
(K G — k) + 3 , etc l34ll . Because the exponent gap is only 
1 /3, it is difficult in simulations to see a clean leading order 
term. Fig.[T9lshows the beating input in Fourier space and in 
physical space. 



Appendix D: The Orr-Sommerfeld problem for Galerkin 
truncation 



Here we study the spectral properties of the Orr- 
Sommerfeld operator governing weak perturbations near the 
time = 1 of the first singularity of the untruncated solution 
for the initial condition uq = sin a;, as introduced in Sec. Hill 
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FIG. 19: (Color online) The beating input term, calculated at and 
used together with the freezing approximation, evaluated for the sin- 
gle mode initial condition Vo(x) = Uo(x) = — sin(a;) (origin at 
preshock) for Kg = 700. (a): imaginary part in Fourier space; (b): 
physical space. 



This operator involves the Fourier transform u i , t k of the solu- 
tion of the untruncated Burgers equation at time t*. Because 
this is an odd pure imaginary function of k, we set u ir> f. = iuk, 
where ilk is a real odd function. With this notation the opera- 
tor becomes the following real (2K G + 1) x (2K G + 1) matrix 



.4 



kk' 



kiik-k 



-K G < k< K a 



-K r , 



<k'< K G . 

(Dl) 

Indices are integers but not restricted to nonnegative values, 
unless otherwise stated. Henceforth, unless otherwise stated, 
all summations are from — K G to K G . 

We studied the spectral properties of A by standard numeri- 
cal techniques using MATLAB® , but we also have some ana- 
lytical results of an algebraic nature. It is appropriate to begin 
with the latter. 

The matrix A is singular, that is it has a vanishing de- 
terminant. This is equivalent to stating that there exists a 
non-vanishing vector (here called a zero mode) ip such that 
5Zfc' Akk'i^k' = for all k. If we can find a ip such that 

k'tp k '=0, for all k, (D2) 



k 1 



then it follows from ( IDU that ip is also a zero mode of A. The 
matrix Uk-k' is skew-symmetric and of odd dimension. By 
an elementary theorem of Jacobi its determinant vanishes and 
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FIG. 20: (Color online) (a): imaginary part of the eigenvalues for 
A'g = 700; (b): ratio of imaginary to real part of the eigenvalues 
for K G = 700; (c): rescaled imaginary parts of the eigenvalues for 
K G = 100, 200, 300, . . . , 1000 showing collapse. 
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FIG. 21: (Color online) Eigenvectors of the Orr-Sommerfeld op- 
erator in Fourier space (origin at the tyger) for K G — 700; real 
parts are shown in red (black) and the imaginary parts in cyan (light 
grey), (a) modes corresponding to an imaginary part of the eigen- 
value which is zero (zero mode, which can be taken purely real), (b) 
halfway between zero and threshold, (c) at threshold, (d) halfway be- 
tween threshold and the largest value, (e) close to the largest value, 
(f) largest value (the inset zooms on the largest wavenumbers). 



thus it has a zero mode. Since the entries Uk-k' are real, the 
zero mode can also be taken real. 

The non-vanishing (complex) eigenvalues come in opposite 
pairs; the associated eigenvectors which are even or odd in k. 
The eigenvalue/eigenvector equation for the Orr-Sommerfeld 
operator A reads 



^kuk-k'tpk' = M'k 



(D3) 



Observe that the operator A dD 1 b is neither Hermitian nor an- 
tihermitian and that, for non-vanishing A, we have tp = 0. 
We now exploit the oddness of iik, to look for even and odd 
eigenvectors. First assume that ipk = ip-k- In ( ID3I ). limit- 
ing ourselves to k > we separate the k' contributions into 
positive and negative ones and obtain, using the oddness of 

k' = K G 

^ k {Hk-k' + Uk+k')ipk' = >4k, k > 0, k' > 0. 



fc'=i 



(D4) 



Now we rescale our eigenvectors by a factor of 1/ Vk 1531 : 

<t>k = (D5) 
so as to rewrite the eigenvalue/eigenvector equation as 

k'=K G 

^/kk 1 (uk-k' +u k +k')(f>k' = A0 fc , fc>0, fc'>0. 

k' = l 

(D6) 

Proceeding similarly under the asumption of an odd eigenvec- 
tor, we obtain instead of ( 




u k +k') 4>k' = Mk, k > 0, k' > 0. 



We now observe that 



(D7) 

and ( ID7l i are two eigen- 
value/eigenvector equations involving two K G x K G "re- 
duced" matrices which are negative transposed of each other. 
Hence their eigenvalues are opposite. 

If all eigenvalues of A and of the two reduced matrices are 
distinct (something for which we have so far only numerical 
evidence), then the even eigenvectors, the odd ones and the 
zero mode exhaust the list of 2A' G + 1 eigenvectors of A. 

The other results on the spectral properties of the Orr- 
Sommerfeld operator are obtained numerically (mostly for the 
case K G = 700) and described now, with occasional soft phe- 
nomenological interpretations of the findings. 

The eigenvalues other than zero are all complex but very 
close to being pure imaginary. Fig. [20] for the case K G = 700, 
shows the imaginary parts of all 1401 eigenvalues: they range 
almost exactly from —K G to K G . Probably, this is related 
to the fact that the Orr-Sommerfeld operator d24l is a modi- 
fied advection operator with an advecting velocity that ranges 
from —1 to +1 and that the k factor which stems from the 
space derivative cannot exceed K G because of the truncation 
operator. 
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The second panel in Fig.[20]shows the ratio of the real to the 
imaginary part, which is quite small, even for the eigenvalues 
close to zero. Actually, the real parts can be neglected alto- 
gether when used in ( f33l >. The explanation may lie in the de- 
composition of the matrix iD6i into the antihermitian matrix 
y/kk'uk-k' an d the Hermitian matrix — ykWuk+k' ■ Remem- 
ber that in this reduced formulation both k and k' are positive. 
In the Hermitian part k, k' and k + k' must be less or equal to 
K G . This prevents k and k' from being simultaneously close 
to K G . But we also found that most of the eigenmodes are 
confined essentially in a relatively narrow boundary layer near 
K G . Hence the Hermitian part cannot contribute much J56ll . 

The third panel in Fig. 1201 shows that by simply rescaling 
by factors cx K G the horizontal and vertical axes of the distri- 
butions of the imaginary parts of the eigenvalues, the curves 
for different K G all nicely collapse on top of each other, for 
large enough K G , suggesting that some limiting distribution 
exists as K G — > oo This we interpret — in highly speculative 
mode — as follows. Because of the confinement near K G of 
most of the eigenmodes, we can approximate the reduced ma- 



trix y/kk'iik-k 1 by K G Uk-k', which is — \K G times the con- 
volution with the Fourier transform of the velocity u+. The 
eigenfunctions of this operator are Dirac measures and the 
corresponding eigenvalues are — \K G times the values the ve- 
locity takes at the supports of these Dirac measures. 

We now turn to the eigenmodes ("Fig. [21~b. It is seen that the 
zero mode (first panel: eigenvalue zero) is nearly constant, 
around 0.0265 with additional small-amplitude even-odd os- 
cillations and some edge effects near ±K G . This structure is 
not surprising: if it was not for truncation effects, (ID2t would 
be the Fourier transform of the equation u+(x)ip(x) = 
whose solution is an arbitrary linear combination of Dirac 
measures at the two zeros of u*, one at the tyger (x = 0) 
and one at the preshock (x = ir). By Fourier transformation, 
these go over into a constant vector and a vector proportional 
to (— l) fe . As we move to higher eigenvalues, we find that 
the eigenmodes are localized at higher and higher wavenum- 
bers. This is unexplained but consistent with the almost purely 
imaginary character of the eigenvalues, as given above. 
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